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coefficient  of  thermal  expansion 
dimensionless  quantity  in  contact  equation 
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Fig.  1.  Idealized  configuration  of  brush.  Outer  surface  R is 
isopotential,  and  the  contact  patch  is  -X  < x < X. 


Z-plane 


Fig.  2.  Axisymmetric  configuration,  with  inner  radius  e, 
and  outer  radius  R. 
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I.  SUMMARY 

Experimental  results  in  the  literature  indicate  that  electric  brushes 
may  deform,  leading  to  point  contact  with  the  slip  ring  at  high  current  * 

levels.  The  present  work  provides  theoretical  support  for  these  observations 
and  a quantative  statement  showing  how  operating  variables  influence  this 
behavior.  The  analysis  is  restricted  to  a two  dimensional  model  of  the 
brush,  with  isotropic  and  constant  material  properties.  The  most  unexpected 
result  is  that,  for  steady  current  flow,  cooling  of  the  exterior  of  the 
brush  increases  the  tendency  to  deform  unfavorably. 


II.  INTRODUCTION 


The  electrical  brush  may,  in  idealized  form,  be  thought  of  as  a 
I wafer  shaped  object  making  contact  with  the  moving  surface  of  a steadily 

turning  cylindrical  body,  or  slip  ring.  The  line  of  contact  would  be 
parallel  with  the  axis  of  the  cylinder,  and  the  sliding  would  be  perpen- 
dicular to  this  line.  Even  when  no  current  passes  through  the  line  of 
contact,  frictional  heating  may  deform  the  brush,  causing  contact  to 
transform  from  a uniformly  loaded  line  to  one  or  more  discrete  patches. 

This  is  the  result  of  the  formation  of  thermal  asperities  which  move 
slowly  across  the  contact  zone.  These  thermal  asperities  are  often 
several  multiples  of  the  height  of  the  initial  roughness  or  waviness  of 
the  brush  face,  and  the  peak  of  each  corresponds  to  a small  region  or 
patch  of  contact  with  the  slip  ring,  which  may  be  highly  stressed  as 

well  as  hot.  Several  studies  have  clarified  the  nature  of  the  transition 

1 2 

to  patch-contact  in  the  absence  of  current  flow,  ’ and  Dow  has  verified 

3 4 

the  theoretical  predictions  experimentally.  Kilaparti  has  investi- 

gated the  role  of  wear  and  other  factors  at  work  in  the  limited  zone  of 
contact . 

Approaching  the  problem  from  the  electrical  side,  Marshall  has 

shown  that  a similar  formation  of  heated  contact  patches  occurs  in 

5,6  7 

brushes  carrying  high  levels  of  current.  McNab  has  reviewed 

the  problem  of  brush  wear  and  has  noted  ambiguities  in  the  literature 
as  to  the  combined  effects  of  current,  load  and  speed.  Undoubtedly 
these  interactions  are  complicated  by  the  formation  of  patch  contact. 

The  present  study  is  intended  to  provide  an  analysis  which  incorporates 
both  frictional  heating  and  current  flow,  and  establishes  how  these  and 
other  factors  interact  to  determine  contact  patch  size,  stress  level  and 
temperature . 
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III.  STATEMENT  OF  THE  PROBLEM 


I 


► 

■ 


I 

A two  dimensional  brush  of  unit  thickness  will  be  assumed,  of 
roughly  the  configuration  shown  in  Fig.  1.  It  is  approximately  in 
the  form  of  half  a circular  disk  and  has  a flat  side  or  edge.  Near 
the  center  of  the  edge  is  the  contact  patch  of  width  2^.  The  outer 
circular  boundary  will  be  taken  as  an  isopotential  line  over  which 
a uniformly  distributed  current  passes  into  the  brush.  Current  exits 
through  the  contact  patch.  Assuming  the  electric  field  is  quasistatic 
it  will  obey  LaPlace's  equation;  and  once  the  potential  distribution 
is  determined,  the  distribution  of  electrical  heating  may  be  found. 

Hie  resultant  temperature  field  can  be  obtained  in  two  steps.  First 
a solution  is  obtained  for  the  condition  that  electrical  heat  only  flows  out 
through  the  contact  patch,  and  that  all  other  boundaries  are  adiabatic. 

Next  a solution  is  found  for  heat  flow  into  the  contact  patch  and  out 
through  the  curved  surface  of  radius  R.  By  superposing  these  two  solutions, 
suitably  scaled,  any  combination  of  heat  flows  through  the  patch  and  outer 
surface  can  be  produced. 

Once  the  temperature  fields  are  specified,  the  thermal  displacement 
of  the  boundary  can  be  calculated.  In  many  instances  a rounded  bulging- 
out  of  the  contact  patch  will  be  predicted,  and  through  application  of 
Hertz's  theory  of  contact  of  cylindrical  bodies  the  conditions  can  be 
found  where  the  contact  patch  would  indeed  be  pressed  smoothly  against 
the  slip  ring  surface  with  a contact-free  gap  to  either  side. 
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IV.  METHOD  OF  ATTACK 


* 

i 


Because  the  thermal  and  electrical  problems  can  easily  be  solved 
in  an  axisymmetric  configuration  (see  Fig.  2) , it  has  been  found  practical 
to  transform  these  solutions  into  a very  close  approximation  of  the  chosen 
configuration,  by  conformal  transformation.  Once  the  temperature  fields 
are  so  generated,  a simple  integral  equation  can  be  applied  to  determine 
the  deflection  of  the  edge  of  the  wafer. 

A Green's  function  for  edge  deflection  can  be  generated  through 
determination  of  the  deflection  distribution  produced  by  a small  patch 
of  elevated  temperature  at  an  arbitrary  position  in  a semi-infinite  plate. 
This  can  be  found  by  adaptation  of  the  equations  for  stress  and  deflection 
in  a heated  axisymmetric  body  in  plane  stress.  In  Fig.  3 is  shown  the 
body  with  the  uniform  temperature  patch  of  radius  e.  Disregarding  for 
the  moment  the  broken-line  construction  also  on  the  figure,  we  note  that 
for  points  external  to  the  patch  the  stresses  and  deflections  are  given  by: 

a = -Q'ETe2/2(x2  + y2) 
o0  = c*ETe2/2(x2  + y2) 

ur  = (1-Kj)  c*Te2/2Vx2  + y2  (1) 

Using  the  method  of  images  as  shown  by  the  broken  lines  in  Fig.  3, 
we  may  superpose  the  fields  produced  by  two  temperature  sources  equidistant 
from  the  x-axis,  and  find  that  on  the  x-axis: 

a = 2cr  cos  20 
y r 

a =o 
xy 

U = 0 (2) 

y 

For  derivation  of  (1)  and  (2)  see  Appendix  A. 
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By  application  of  the  integral  equation  derived  in  Ref.  (9)  for  edge 
displacements  of  a plate  having  an  edge  load  o , we  find  the  normal 
displacement  of  the  edge,  6,  to  be: 


6 = 4:  l CTy(5)Xn|S-x|dS 


(3) 


where  5 is  a dummy  variable.  Integrating  and  letting  ne  = dA^ , one  finds 
that  6^  is  the  deflection  of  the  plate  edge  at  x = 0,  caused  by  a temperature 
patch  T^dA^ , at  the  position  x^  , y^. 

dA. 


40'T,  fTTy  , 

6i  i~T  + lx1lA»ly1/*1l} 


xi  + yi 


When  there  is  a continuous  distribution  of  temperature  in  the  body  the 
influences  may  be  summed  to  give,  see  Appendix  B 
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Ua  f Tny  A 
— u L 2 + 

TT  A 


•^n  |y/x  | 


TdA 


, 2 2 

(x  + y ) 


(4) 


As  it  stands  this  equation  gives  the  displacement  at  the  origin  of  the 
coordinate  system,  but  a simple  geometric  transformation  will  cause  it 
to  apply  at  x = 5. 


‘w-^f  Ji>+  irtWrilil 

TT  A 


TdA 


y“+(x-£)2 


(5) 


it  is  of  practical  interest  that  the  second  term  in  the  brackets  has  only 
a small  influence  on  relative  displacements  in  the  contact  patch. 


V.  TEMPERATURE  SOLUTIONS  IN  THE  AXISYMMETRIC  PLATE 

For  the  case  of  simple  conduction  through  the  cylindrical  surface  of 
radius  e,  see  Fig  2,  LaPlace's  equation  for  temperature  will  be  satisfied 
by 

T - T(r)  = (Q/KTT)4n(r/e)  (6) 

c b 

where  Q,  is  the  heat  flow  moving  radially  outward,  T is  the  temperature 
b c 

of  the  inner  cylindrical  surface  and  T(r)  is  the  temperature  at  arbitrary 
r interior  to  R. 

Turning  now  to  the  question  of  electric  current  flow,  one  finds  by 
analogy 

e - e(r)  = --jr  4n(r/e)  (7) 

C TT 


Noting  that  electric  heat  generation  will  be  given  by 


1 ,de. 2 

« = p'(d7) 


(8) 


it  follows  that  the  total  electrical  heating  between  r and  R is  given  by 


Q(r)  - ^n(^) 


(9) 


irrespective  of  the  direction  of  current  flow.  The  maximum  magnitude  of 
Q would  be 

2 

Q,«V  MR/*)  d°) 

'e*  m 


and  would  be  realized  as  heat  flow  through  the  inner  surface,  if  all 
heat  flow  were  blocked  at  the  boundary,  R.  Under  the  same  boundary 
condition,  we  may  write  at  any  r, 


Q(0 


_ I JL 

TT 


R 

r 


6 
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and  it  follows  that 


T(r) 


T(R) 


[in(r/R) ]2 
2Kn  in(R/e) 


This  may  be  rewritten  as 


T(r)  ‘ T(R)  = {in<R/£>  + 2Xn(e/r)  + 

xn(R/e) 


(ID 


(12) 


i 


VI.  CALCULATION  OF  SURFACE  DEFLECTION 

As  stated  above,  solutions  in  the  W-plane  were  to  be  obtained  from 
the  axlsynunetrlc  Z-plane  solutions  by  conformal  transformation.  To  this 
e...  the  Zhukovsky  transform  may  be  applied. 

W = Z + 1/Z  (13) 


This  may  be  rewritten  as 


Z = [ W + V w2  - 4 1 


/ 2 


(14) 


Hence  for  any  W,  expressed  as  a complex  number  x + iy,  there  will  be  a Z, 

iV 

which  can  be  expressed  as  re  . Knowing  r,  one  may  use  Eq . (6)  or  (12) 
as  appropriate,  to  obtain  the  temperature  at  the  point  and  at  the  corre- 
sponding W-point.  The  use  of  the  transform  on  the  electrically  heated 

field  is  permissible  because  the  quantity 

1 ~ i de\2  de^' 

q = i LAaf)  +vs^J  (15) 


transforms  such  that 


q(Z)  = q(W) |4|!' 


(16) 


See  Ref.  (10).  Once  values  of  T are  known  at  points  in  the  W-plane  it 
is  possible  to  evaluate  the  integral  of  Eq . (5)  at  various  values  of 
Closed  form  solution  was  possible  only  for  the  case  of  simple  heat  transfer 
(see  Eq.  (6)).  Consequently,  the  space  was  subdivided  into  a net  of  points 
and  the  influences  of  each  were  summed.  The  simple  heat  flow  solution 
provided  a check  on  the  accuracy  of  the  summation  procedure.  Results  are 
reported  here  as  the  difference  in  displacement  between  contact  center 
and  edge  or 

A 

6 = 6(0)  - 6(X)  (17) 
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For  the  case  of  simple  heat  flow  radially  outward,  the  exact  solution  is: 


k 


*«Qh 

Sth  * "ttk  <°-5708> 


(18) 


By  numerical  integration  (for  R/e  ■ 1000) 

/oQ, 

6 “ 

See  Appendix  C. 


TTK 


(0.5715) 


(19) 


The  close  agreement  gives  confidence  in  the  numerical  procedure,  so  it 
may  be  applied  to  the  case  of  electrical  heating  with  all  of  the  heat 
passing  outward  through  the  contact  patch,  drawing  upon  Eq.  (12)  for 
the  temperatures,-  and  noting  that  any  component  of  the  temperature  distri- 
bution that  is  constant  throughout  the  field  will  not  contribute  to  defor- 
mation of  the  boundary,  therefore  it  will  not  contribute  to  $.  The  result 
of  this  numerical  integration  is 


l0fQe  r 


0.7636  1 


TTK  i_‘0,5708  + I^rT®)] 


(20) 


This  result  is  interesting  in  two  respects.  First  it  represents 
an  indentation  rather  than  a protrusion  of  the  surface,  and  second  it  is 
dependent  on  the  size  of  the  brush,  R.  In  the  first  we  note  simply  that 
for  heat  to  be  removed  through  the  contact  patch  the  lowest  temperature 
must  be  at  the  contact  surface,  hence  a deficit  of  temperature  there  and 
the  resultant  indentation. 

The  second  effect  is  the  result  of  the  fact  that  the  more  material 
through  which  current  is  passed  the  more  heat  generated. 

The  combined  displacement  for  electrical  heating  and  brush  cooling 
may  be  obtained  by  superposition  of 


or 


/ofQ 


TTK 


(21) 


Since  represents  heat  passing  through  the  outer  surface  of  the  brush 
it  may  be  thought  of  as  a measure  of  cooling.  If  the  brush  were  insulated 
would  be  zero.  With  some  forced  cooling  /Qe  might  be  any  value  up  to 


I 


or  exceeding  unity. 
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VII.  COMBINED  EFFECTS  OF  FRICTIONAL  AND  ELECTRICAL  HEATING 


Frictional  heating  of  the  contact  patch  is  given  by 


Qf  = M.VP 


where  V is  sliding  speed,  P is  contact  load  and  M-  is  friction  coefficient. 
Recalling  that 


Qe  = (I2p  /rr)Xn(2R/i 


(22) 


one  may  write 


‘total 


= + Q, 


(23) 


We  now  note  that  with  isolated  slip  ring  and  cooled  brush  could  have 
as  its  upper  limit  Qtotaj_*  Hence  the  maximum  value  of  6,  given  by  Eq . (22) 
would  be 


6 

max 


^L0.57°8(Qtot-Qe 


) + 


0,7636 
^n(R/e)Qe  - 


(24) 


For  convenience  call  the  quantity  in  brackets  Y,  (in  Eq . (24))  making 

6 = (aX/nK)(y)  (25) 


Let  us  now  draw  upon  the  well  known  relationship  for  Hertzian  contact 


6 = P/1.72E 


(26) 


A 

where  6 is  the  distance  the  center  of  the  contact  is  indented  relative 
to  the  edges  of  the  contact.  It  follows  that  if  the  thermal  bump  is 
pressed  flat,  to  assure  full  contact  over  the  region  2X, 


P 

1.72E 


(27) 


10 


11 


Here  it  is  assumed  that  the  deflection  of  the  slip  ring  is  negligible 
relative  to  that  of  the  carbon,  therefore  only  the  Young's  modulus  E 
of  the  carbon  brush  need  be  considered. 

In  Eq . (27)  it  is  seen  that  & increases  with  brush  load  P.  Increasing  ) 

of  cooling  of  the  brush  causes  Y to  increase  and  causes  ^ to  be  smaller  for 
a given  load,  an  adverse  condition.  Although  the  derivation  is  not  valid 
for  large  & relative  to  brush  size  R,  nevertheless  the  condition  = R 
provides  a criterion  as  to  when  thermal  effects  are  beginning  to  be  important, 
and  for  X«  R the  equation  should  serve  to  predict  contact  patch  size. 

To  obtain  contact  temperature  one  need  only  to  refer  to  Eq.  (5,12) 
and  superpose  the  magnitudes  of  T(r  = e) , thus  giving  the  elevation  of 
contact  temperature  above  brush  temperature  T(R) . To  obtain  contact  stress 
note  that  its  mean  value  would  be  P/2^. 
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APPENDIX  A 


DERIVATION  OF  STRESSES  AND  DEFLECTIONS  FOR  POINTS  EXTERNAL  TO  THERMAL  PATCH 


From  reference  1,  one  can  find  the  expressions  for  a thin  circular 


disk  with  the  edge  r ■ b free  from  stress: 


1 rb  1 

f A m 1 


- crECg^  Trdr  - 1 Trdr) 

o 


2 J 

r o 


,b 


i rr 


<*0  = QfE(— ^ i Trdr  + — | Trdr  - T) 


.2  J 
b o 


2 - 
r o 


(A-l) 


1 rr  r 'b 

u = (1+v)q  — j Trdr.  + (l-v)or  — J Trdr 
o b ” o 


For  a body  with  a uniform  temperature  patch  of  radius  e,  the  integral  in 
the  brackets  is 


rr  ,-b  Te2 

I Trdr  = Trdr  + Trdr  = — 

v J i Z 

o o e 

•r  'r  Te2 

Trdr  = Trdr  + Trdr  = — r— 

%•  %i  %)  z 

O O e 


(A-2) 


Where  the  second  integral  is  zero  when  the  temperature  is  zero  outside  of 

2 

patch.  And  note  that  ne  = dA. 

For  b»e. 


_ W1  1 . TdA 

CTr  = “E(b  - ~ 2 IFT 
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aETdA 

2 2 
2tt(x  +y  ) 


1 TdA  . 1 TdA  T)  _ QfETdA 


0 3 ffE(,2  2n  + 2 2n 
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u = (l+v)QrTdA/2TT.  x2+y2 
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Using  the  method  of  images  as  shown  by  the  broken  lines  in  Fig.  3, 
we  may  superpose  the  fields  produced  by  two  temperature  sources  equi- 
distant from  the  x-axis,  and  note  that  by  Mohr's  circle  one  can  find  that 

CTy  = 2ar  cos  20  (A-3) 

CT  = 0 

xy 

u = 0 

y 


APPENDIX  B 


DEFLECTION  OF  THE  PLATE  EDGE 


From  (3)  one  finds  that  the  deflection  of  the  plate  edge  at  x = 0 is 

from  (2) 

a (§)  = 2o  cos  20 
yi  r 


(B-l) 


= 2<7r(l-2  cos  0) 


, 2f  - aETidAi  Ti 
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p2  2 ‘ .-2  2 2-j 

5 +y.  (?  +yt) 


(B-2) 


Substitute  (B-2)  into  (B-l)  and  note  that  deflection  is  symnetric  to  y-axis. 
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(B-3) 


■x  in(x  -5)  ^ in(C-x.) 
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One  can  integrate  each  term  in  (B-4)  by  parts  as  following: 

f‘  ta„-i  i.  |*i  + x f* 
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After  substitution  one  can  arrive. 
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When  there  is  a continuous  distribution  of  temperature  in  the  body,  the 
influence  may  be  integrated  over  the  whole  body  to  give 


IVA.^r?+uu„|*|;^ 


TdA 
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(B-5) 


APPENDIX  C 


EXACT  SOLUTION 


For  the  case  of  simple  conduction  the  heat  flow  through  the  contact  line 
segment  in  W-plane  is 


I dZ| 

qw  = VdW1 

where  qw  and  q ^ are  heat  flow  per  unit  length  in  W-plane  and  Z-plane. 
If  we  let  e equal  to  unity  in  Z-plane,  then 

q = 2k  = fib 
ire  tt 

From  (14)  one  can  get 


(C-l) 


(C-2) 


m = 

dW 


2 ;1-(|)2 


(C-3) 


the  semi  unit  circle  transforms  into  a line  segment  in  W-plane.  Also, 
one  can  find  that11  the  curvature  in  W-plane  is 

dV  -°% 


dx 


,2 


(C-4) 


here  for  qw  heat  flow  moving  outward,  the  curvature  is  negative. 


let  x'  = £ x'U  1 and 


-^z.dZ, 

dx2  K ' K 'dW1 

-2ttQk 


KTT  VI 
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(C-5) 


integrate  (C-5)  twice  and  note  dimensional  deflection  6=6'^  one  can  get 


6(x')  » 

here  c is  a Constance. 


ttk 
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1 / 2 
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appendix  d 


NUMERICAL  SOLUTION 


From  (5)  and  (6)  we  have 


dA 


note  that 


n“  A 
2 

dW. 


2 2 

y +<*-5> 


% - ".if! 


4 2 

r - 2r  cos  20+1 


drd0 


(D-l) 


x = (r  + i)  cos  0 
r 

y = (r  - i)  sin  0 

one  can  integrate  from  r=l  to  R 0<=O  to  n numerically,  and  note  that 
dimensional  6 = - 6'.  The  final  result  is 

5 = 6(0)  - 6(2)  = (0.5715)  (D-2) 
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COMPUTER  LISTING 


PROGRAM  l‘(T  (O'lTPUT) 

c numcpical  if.rror.'AT low  or  oEFUCTION 

C A = A'lOl.E  FROM  n TO  1,U(|6 

C RADIUS  lMcr?EMF,MT  = G\ 

C A.'v'GULAF  lMCnrMrMT  = G ’ 1 • 7 0 >'*. 

n = 1 OV'.  - 
G = o . " 0 ? 

P = 0,1 
1 o = O.o 

ON  = 0.0*1  IS 
SP  = 1.0  ♦ G A / P.O 
? A = 0 * 1.S73H  / p.r 

3 A?  = A ♦ ;\ 

SI  = vlhi  (A) 

Cl  = COS  (A) 

C ' = ros  (A?) 

T=  >>*SR*SR*r  P.  1 . ->)  / RP**3.»i 

111  = <SP  4 l.-‘/SK  1 « rl 

I V = <$n  -!.'•/  r>R  ) * '»1 

ll?  s :||  - B 

tl  - APS  (UP) 

S = U / V 

Q=(3.  l4U>*V/P-ll*ALOCi(S>  ) / (IJ»IJ  ♦ V«V) 
F 1 = P / _M 
F = AlOG  (I'M 

. pu  = r » T » .j 

ON!  = pM  * (-  « l • r> 7 0 H * ON 
n = o ♦ l>i 
A = A ♦ d * 1,070* 

IF  ( A .IF.  3.1410  ) nr,  Tf  3 
TF  ( ? ."".  O.i  ) 00  to  /, 

PRli.lT  10>i  , Ok,  (,A 

4 SP  = CO  ♦ r-A  / ?,'■ 
r,A  = !.?  u nr. 

— - • - SR  = cp  ♦ r-A  / p. a 

IF  (SP  .IF.  P ) 0.1  TP  ■’ 

TF  (?  .Mr.  0,0  ) ro  TP 
DC  = 3 

?'  x = ( P - DC  ) <>  2.0  / 3.1410^3  **  ?,n 

PRINT  ?(lO  » P,  P • X . n 

. P = P ♦ .?  . P 

IF  ( P .L*.  2.0  ) on  TP  1 
p*r.  FORMAT  (2(FX,F4.D  # 2 (S<  *f  20  ,/>)  1 
_.V.t  • FORf  j AT  (2(1  OX. Flu. 4)  ) 

STOP 
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